Strain-induced ferroelectricity in CaTiOa from first principles 
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First principles calculations are used to investigate the effects of epitaxial strain on the structure 
of the perovskite oxide CaTiOs, with particular focus on the stabilization of a ferroelectric phase 
related to a polar instability hidden in the orthorhombic equilibrium bulk Pbnm structure but found 
in previous first-principles studies of the ideal cubic perovskite high-symmetry reference structure. 
At 1.5% strain, we find an epitaxial orientation transition between the ah-ePbnm phase, favored 
for compressive strains, and the c-ePbnm phase. For larger tensile strains, a polar instability 
develops in the c-ePbnm phase and an epitaxial-strain-induced ferroelectric phase is obtained with 
polarization along a <110> direction with respect to the primitive perovskite lattice vectors of the 
square substrate. 
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With the recent dramatic advances in the synthesis of 
coherent epitaxial films of complex oxides [ij , it is possi- 
ble to maintain extremely high strains in thin film materi- 
als, often as much as 2-3%. In some materials, such high 
strains can drive the system through a structural phase 
boundary to a novel phase with structure and properties 
distinct from those of the bulk equilibrium phase. First 
principles calculations can provide quantitative predic- 
tions about these novel phases and phase boundaries; in 
particular, such phases can be identified from examina- 
tion of the lattice instabilities of the high-symmetry ref- 
erence structure of the bulk phase [2] , and their structure 
and properties predicted. 

A prototypical example is strain-induced ferroelectric- 
ity in SrTiOs. Starting from the bulk equilibrium para- 
electric phase, strain-polarization coupling results in a 
polar instability beyond critical values of both tensile 
and compressive (001) epitaxial strain. This behavior 
was predicted by Landau theory 0] and further ana- 
lyzed in first principles investigations 0, . Experimen- 
tal observation of epitaxial-strain induced ferroelectricity 
in SrTiOa 0, Q demonstrates this strain-induced ferro- 
electricity both for compressive and tensile strain. 

CaTiOa presents a greater challenge for the observa- 
tion of strain-induced ferroelectricity. A relatively large 
temperature-dependent dielectric response has led to the 
characterization of CaTiOa, like SrTiOa, as an incipi- 
ent ferroelectric However, unlike SrTiOa, CaTiOs 
has a strong tendency to oxygen-octahedron rotations, 
which tend to suppress polar lattice instabilities [1, Q. 
Indeed, the bulk Pbnm structure [lo| . a structure type 
that includes a large number of other perovskite oxides, 
is obtained as a distortion of the high-symmetry cubic 
perovskite reference structure by freezing in components 
of the M+ and R| oxygen octahedron modes (notation 
from Ref. ITU ) involving rotation around [001] and tilt- 
ing around [110], respectively, with additional changes 
in the lattice constants a, b and c and Ca displacements 
preserving the space group symmetry. All other known 
phases are also nonpolar: these include a paraelectric cu- 



bic structure at high temperatures and two intermediate 
phases (one tetragonal and one orthorhombic) [l^- Less 
is known about high-pressure phases, though a transi- 
tion to a nonpolar orthorhombic Cmcm structure and, 
at higher pressures, to a nonpolar post-perovskite struc- 
ture like that of MgSiOs have been proposed based on 
first-principles results ^ISj. 

First-principles calculations of the full phonon disper- 
sion relation of the cubic perovskite structure of CaTiOs 
show that in addition to the expected instability of M;^ 
and R4 , the cubic perovskite structure does have a 
strongly unstable polar F4 mode [l3, [3, Freez- 
ing in of this polar mode, which involves displacements 
of the Ca and Ti ions relative to the oxygen octahedron 
network, would yield a ferroelectric phase with nonzero 
spontaneous polarization. The fact that it does not con- 
tribute to the observed bulk phases suggests that it is 
inhibited by the oxygen octahedron rotations. If the ro- 
tations are artificially suppressed, as is possible in a first- 
principles calculation, the FJ mode dominates and the 
resulting ferroelectric PAmm phase is found to have a 
very large polarization 
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To stabilize a ferroelectric phase of CaTiOs under con- 
ditions realizable in the laboratory, it is necessary to 
change the balance of the competition between the oc- 
tahedral rotations and the polar instability in favor of 
the latter. As in SrTiOa, the polar instability in CaTiOs 
could be strengthened through the well-established sen- 
sitivity of the polar mode to strain in the titanates 17| , 
specifically by tuning epitaxial strain 3, 3- However, 
the oxygen-octahedron distortions inevitable in CaTiOs 
will couple to the epitaxial strain and polar instability, 
leading to a richer set of possibilities than in SrTiOs. 

In this paper, we study the use of epitaxial strain to 
stabilize a ferroelectric phase of CaTiOs. We focus on the 
orthorhombic bulk ground state, and investigate whether 
epitaxial strain can induce a polar instability, analogous 
to the behavior of SrTiOs. The relatively low symme- 
try of the Pbnm structure requires careful attention in 
imposing the epitaxial constraints, and introduces new 
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TABLE I: The Wyckoff parameters of the Phnm ground state 
and the c-ePbnm structure at 3% tensile strain. 



Phnm 



c-ePhnm 



Ca (4c) X = 0.510, y = -0.047 x = 0.509, y = -0.042 
Ti (4a) 

O (4c) x = 0.081, y = 0.021 a; = 0.085, y = 0.022 

X = 0.207, y = 0.292, x = 0.218, y = 0.283, 

z = -0.043 2 = -0.044 
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features into the strain-induced ferroelectric state. 

We performed density functional theory total-energy 
calculations within the LDA approximation using the 



VASP code 
tials 
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22] with its supphed PAW poten- 
23l , |24| . The plane- wave energy cutoff was 680 cV 
and the k-point grid was 8x8x6. 

For the study of the effects of epitaxial strain, we 
carried out "strained bulk" calculations, in which total- 
energy calculations are performed for the periodic crys- 
tal with appropriate epitaxial constraints imposed on the 
lattice parameters. In some cases, these are constraints 
which cannot be automatically imposed within the avail- 
able VASP relaxation algorithms. In such cases, we de- 
veloped an elastic energy expansion around the lowest 
energy Pbnm structure that satisfies the epitaxial con- 
straint by fitting to the energies of structures with small 
changes in strain (the latter structures not necessarily 
satisfying the epitaxial constraints). This elastic energy 
is then minimized with respect to strain subject to the 
epitaxial strain constraint; the resulting lattice parame- 
ters are then fixed in a total-energy calculation in which 
the internal structural parameters are relaxed. All struc- 
tures were relaxed until the forces on the atoms were less 
than 2.5 meV/A. 

For selected structures, we compute the stability 
against zone-center modes by performing frozen phonon 
calculations in which single atoms are displaced by ap- 
proximately 0.01 A. From finite differences of the result- 
ing forces, the force constant matrices are determined and 
subsequently diagonalized to obtain eigenfrequencies and 
eigenvectors. 

Our results for the structure of the bulk orthorhombic 
ground state are given in Table |T1 The structure has an 
energy 410 meV/f.u. lower than the ideal cubic perovskite 
structure. Consistent with previous first-principles cal- 
culations [l^, [3] , we find good agreement between the 
computed structure and experiment, taking into account 
that in the local density approximation lattice constants 
typically tend to be underestimated by about one per- 
cent. 

Next, we investigate the effects of epitaxial strain on 
the Pbnm phase. As in Ref. 0, we designate the strained 
phases as ePbnm, where the prefix e denotes "epitax- 
ial." We consider epitaxial strain on a square lattice 




FIG. 1: The two distinct relative orientations of the lat- 
tice vectors and the primitive perovskite substrate matching 
planes in the Pbnm structure are shown for (a) the ab-ePbnm 
phase and (b) the c-ePhnm phase. 



substrate, corresponding to a (001) perovskite surface. 
To allow direct comparison with experiment despite the 
lattice constant underestimate discussed above, we de- 
fine epitaxial strain relative to ao = 3.77 A, which is 
the cube root of the computed volume per formula unit 
of the relaxed Pbnm structure. In the Pbnm structure, 
there are two symmetry-inequivalent primitive perovskite 
(001) planes, as shown in Fig. [1] Thus, there are two dis- 
tinct orientations for an epitaxial film: the first, with tc 
in the matching plane and ta and t^, out of the plane [ab- 
ePbnm, Figure [1] (a)), and the second, with normal to 
the matching plane (c-ePbnm, Fig. [Hb)). 

We compute the total energies for these two orienta- 
tions for epitaxial strains ranging from —1.5% to 4% f25| . 
For c-ePbnm, the c lattice parameter and internal struc- 
tural parameters are relaxed at each strain, maintaining 
the Pbnm symmetry; the a and b lattice parameters are 
fixed by the constraint. The epitaxial constraint allows 
for tc not to be normal to the matching a6-plane and 
tilting tc could lower the energy. However, an elastic 
analysis for the —1.5% and 4% cases showed that tilting 
tc does not lower the energy and we assumed this to be 
true for the intermediate strains as well. ab-ePbnm has 
lower symmetry than c-ePbnm; that is, distinguishing 
one of the two (110) planes removes space group symme- 
tries, resulting in a space group P2i/to; the constraint 
fixes \tc\ and | — ta | , as well as the condition tc - {t\j — ta) 
— 0. To optimize the lattice parameters for this case, we 
use the elastic energy expansion method described in the 
previous section. 

The results are shown in Fig. [2l ab-ePbnm is favorable 
for compressive strains and c-ePbnm is favored with in- 
creasing tensile strains. Within this subspace of nonpolar 
structures, there is an epitaxial orientation transition at 
1.5%. 

Next, we turn to the stability of the ePbnm phases 
against symmetry-breaking distortions, with special at- 
tention to polar phonons. Previous computation of the 
phonon frequencies for the bulk equilibrium structure 
showed three low- frequency polar phonons at 94 cm~^. 
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FIG. 2: Total energy per five-atom formula unit for vari- 
ous epitaxially constrained structures as a function of misfit 
strain. At each strain, the energy of the c-ePbnm structure is 
taken as the zero of energy. The connecting lines are a guide 
to the eye. 



cm ^ and 89 cm ^, with induced polarizations along 



a, h and c, respectively |14| . These phonons are expected 
to be sensitive to changes in strain, based on known 
polarization-strain coupling in calcium titanate 27 1. 

We first consider c-Pbnm with 4% tensile strain, with 
computed structural parameters reported in Table 1. A 
zone-center frozen phonon computation for this struc- 
ture showed four unstable phonons, the lowest two, at 
213i cm"^ and 209i cm^^, being polar and generating 
structures with space groups Pmn2i (polarization along 
a) and Pmc2i (polarization along 6), respectively. In 
both cases, the orientation of polarization with respect 
to the primitive perovskite axes is along the <110> direc- 
tions. The energies of the structures for these two space 
groups, optimized under the epitaxial strain constraint, 
are 35 meV/f.u. and 28 meV/f.u. below c-Pbnm, respec- 
tively, with polarizations 0.46 C/m^ and 0.45 C/m^ com- 
puted using Born effective charges and atomic displace- 
ments [28|. Thus, at 4% tensile strain, we predict strain- 
induced ferroelectricity in CaTiOa. 

For the full range of strains, the unconstrained inter- 
nal structural parameters were optimized within these 
two polar space groups. Only in the 4% case was \tc\ 
re-optimized for the polar structures; this only had a 
marginal effect on the energy and the polarization com- 
pared to when the relaxed value for the c-ePbnm struc- 
ture was used, presumably because the polarizations are 
in the a6-plane. At compressive strain, the nonpolar c- 
ePbnm structure is stable against polar distortions, and 
no ferroelectricity is observed. For tensile strain, the fer- 
roelectric instabilities first appear at 2% strain, and the 
energy gain and polarization of the optimized ferroelec- 
tric phases grow with increasing strain. 

Let us now consider ab-ePbnm with 4% tensile strain. 
We looked for polar instabilities in this structure by dis- 
placing atoms in such a way that the resulting, lowered 



symmetry allowed for nonzero polarization, and then re- 
laxing the internal structural parameters while keeping 
the lattice parameters fixed at their ab-ePbnm values. 
This procedure revealed two structures with polarizations 
in the matching plane of 0.33 C/m^ along c (space group 
P2i) and 0.40 C/m^ along the ab diagonal (space group 
Pm), respectively. (The polarizations were obtained us- 
ing Born effective charges and atomic positions [28j] .) The 
former is the lowest in energy but still above Pmn2i and 
Pmc2i, see Fig. [21 This procedure was also carried out 
for epitaxial strains of 3% and —1.5%. No polar phase 
was found in the latter case. 

The mechanism of strain-induced ferroelectricity in the 
Pbnm phase is closely related to that for SrTiOa, which 
similarly has a high strain sensitivity of the low-frequency 
polar mode. However, the effect in SrTiOa is equally 
strong for compressive as for tensile strain, as elongation 
of the unit cell produced by compressive strain desta- 
bilizes the polar mode as effectively as elongation in the 
in-plane direction for tensile strain. This is a direct conse- 
quence of the fact that SrTiOs is cubic in the paraelectric 
phase; the cubic-tetragonal transition being at low tem- 
peratures and the rotational distortion not being strong 
enough to inhibit the strain-enhanced polar instability. 

In CaTiOa, in contrast, the rotational instabilities are 
much stronger and the resulting distortions are much 
larger. No ferroelectric ePbnm phase was found for com- 
pressive strain, despite elongation of the unit cell along 
the direction normal to the surface; this presumably is 
due to inhibition by the pattern of octahedral rotations. 
This highlights the idea that in a nonpolar low-symmetry 
phase, unlike in a cubic phase, the relationship between 
the crystal axes and the epitaxial constraints is very im- 
portant, different choices yielding quite distinct struc- 
tures and coupling to potential instabilities. 

To investigate the relative importance of enhancing the 
polar instability compared to suppressing the rotational 
instabilities, we analyzed the structural parameters of 
the c-ePbnm phase as a function of epitaxial strain. The 
amplitudes of the and modes change surprisingly 
little in the range of strains reported here, suggesting that 
the dominant mechanism of the strain-induced ferroelec- 
tricity is the strain enhancement of the polar instability. 

Epitaxial strain in the Pbnm is not the only possible 
avenue to a ferroelectric phase. Another promising ap- 
proach is to replace the rotational pattern in the Pbnm 
structure with a different pattern which allows the gain 
associated with the octahedral rotation instability but 
which is less inhibitory to the polar modes; this will be 
discussed in a separate publication [2^ . 

In comparing these predictions with experiments on 
epitaxial strained CaTiOa, it is important to keep in 
mind that our approach considers only the effect of strain 
on the ground state structure and properties. In a real 
thin film, especially the ultrathin films needed to sustain 
very high strains, other factors can affect the observed 
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phase, including temperature, the atomic arrangements 
at the substrate and the film-substrate interface, relax- 
ation, reconstruction, and adsorption at the free surface, 
and defects and impurities in the film itself. However, the 
tendency to ferroelectricity with increasing tensile strain 
is clear in our results, and to the extent that these other 
factors do not act dominantly against it, we expect fer- 
roelectricity in CaTiOs to be observed at experimentally 
accessible strains. Investigations are currently in progress 

In summary, we have performed first-principles cal- 
culations for epitaxially strained structures of CaTiOa. 
At 1.5% strain, we find an epitaxial orientation transi- 
tion between the ab-ePbnm phase, favored for compres- 
sive strains, and the c-ePbnm phase, favored for ten- 
sile strains. For sufficiently large tensile strains, a po- 
lar instability, which was hidden in the equilibrium bulk 
structure, develops in both phases. The epitaxial-strain- 
induced ferroelectric phase lowest in energy originates 
in the c-ePbnm phase and has the polarization along a 
<110> direction with respect to the primitive perovskite 
lattice vectors of the square substrate. 



ACKNOWLEDGMENTS 

We would like to thank M. H. Cohen, V. Gopalan, D. 
R. Hamann, D. G. Schlom, and D. Vanderbilt for useful 
discussions. This work was supported by NSF MRSEC 
DMR-0820404 and ONR N00014-09-1-0302. K. M. R. 
would also like to thank the Aspen Center for Physics, 
where part of this work was carried out. 



[1] D. G. Schlom, L. Q. Chen, C. B. Eom, K. M. Rabe, S. 
K. Streiffer and J. M. Triscone, Ann. Rev. Mat. Res. 37, 
589 (2007). 

[2] K. M. Rabe, in Computer Simulation Studies in 
Condensed-Matter Physics XVI 95, Eds. D. P. Landau, 
S. P. Lewis, and H.-B. Schuttler (Springer- Verlag, Berlin, 
2004). 

[3] N. A. Pertsev, A. K. Tagantsev, and N. Setter, Phys. 

Rev. B 61, R825 (2000). 
[4] A. Antons, J. B. Neaton, K. M. Rabe, and D. Vanderbilt, 

Phys. Rev. B 71, 024102 (2005). 
[5] A. Vasudevarao, A. Kumar, L. Tian, J. H. Haeni, Y. L. 

Li, C.-J. Ekiund, Q. X. Jia, R. Uecker, P. Reiche, K. M. 



Rabe, L. Q. Chen, D. G. Schlom, and V. Gopalan, Phys. 

Rev. Lett. 97, 257602 (2006). 
[6] J. H. Haeni, P. Irvin, W. Chang, R. Uecker, P. Reiche, Y. 

L. Li, S. Choudhury, W. Tian, M. E. Hawley, B. Craigo, 

A. K. Tagantsev, X. Q. Pan, S. K. Streiffer, L. Q. Chen, 

S. W. Kirchoefer, J. Levy, and D. G. Schlom, Nature 

(London) 430, 758 (2004). 
[7] V. V. Lemanov, A. V. Sotnikov, E. P. Smirnova, M. Wei- 

hnacht and R. Kunze, Sol. St. Comm. 110, 611 (1999). 
[8] W. Zhong and D. Vanderbilt, Phys. Lett. 74, 2587 

(1995) . 

[9] N. Sai and D. Vanderbilt, Phys. Rev. B 62, 13942 (2000). 
[10] K. Parlinski, Y. Kawazoe, and Y. Waseda, J. Chem. 

Phys. 114, 2395 (2001). 
[11] H. T. Stokes, E. H. Kisi, D. M. Hatch, and C. J. Howard, 

Acta Cryst. B58, 934 (2002). 
[12] V. Zelezny, E. Cockayne, J. Petzelt, M. F. Limonov, D. 

E. Usvyat, V. V. Lemanov, and A. A. Volkov, Phys. Rev. 

B 55, 224303 (2002). 
[13] X. A. Wu, et al., Sol. St. Comm. 136, 416 (2005). 
[14] E. Cockayne and B. P. Burton, Phys. Rev. B 62, 3735 

(2000). 

[15] R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 49, 
5828 (1994). 

[16] S. M. Nakhmanson, K. M. Rabe, and D. Vanderbilt, 

Phys. Rev. B 73, 060101 (R) (2006). 
[17] R. E. Cohen, Nature 358, 136 (1992). 
[18] S. M. Nakhmanson, K. M. Rabe, and D. Vanderbilt, 

Appl. Phys. Lett. 87, 102906 (2005). 
[19] A. Zayak, S. Nakhmanson, and K.M. Rabe, 2006 

Workshop on Fundamental Physics of Ferroelectrics, 

|httpI7/ www.gl.ciw.ed u/^cohen/meetings/ferro2006/] 

ProgramFerro06 / final.pdf 
[20] Versions 4.6.28 and 4.6.26 of VASP were used. 
[21] G. Kresse and J. Hafner, Phys. Rev. B 47, R558 (1993). 
[22] G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 

(1996) . 

[23] P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

[24] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

[25] The strain values are -1.51%, -1.01%, 1.01%, 2.02%, 
3.03%, and 4.05%. 

[26] A. T. Zayak, X. Huang, J. B. Neaton, and K. M. Rabe, 
Phys. Rev. B 74, 094104 (2006). 

[27] O. Dieguez, K. M. Rabe and D. Vanderbilt, Phys. Rev. 
B 72, 144101 (2005). 

[28] W. Zhong, R.D. King-Smith, and D. Vanderbilt, Phys. 
Rev. Lett. 72, 3618 (1994). 

[29] C.-J. Ekiund, C. J. Fennie and K. M. Rabe, in prepara- 
tion. 

[30] E. Vlahos, C. Brooks, C.-J. Ekiund, K. M. Rabe, V. 
Gopalan and D. G. Schlom, unpublished. 



